import pandas as pd
from osgeo import gdal
import geopandas as geo
import matplotlib.pyplot as plt

# pandas显示所有行
pd.set_option('display.max_rows', None)
# 设置pandas显示的列数
pd.set_option('display.max_columns', 20)
# 让dataframe打印不换行
pd.set_option('display.width', 5000)

# 设置画板
fig = plt.figure(1, (16, 9), dpi=120)
ax = plt.subplot(111)
plt.sca(ax)

# read newYork map
map = geo.read_file('../data/map/NewYork.shp')
# draw newyork map to artboard
map.plot(ax=ax, edgecolor=(0, 0, 0, 1), facecolor=(0, 0, 0, 0), linewidths=0.5)

# read the grid containing only site
grid = geo.read_file('../data/grid/grid_merge.shp')
# print(grid.head())
# draw grid to artboard
grid.plot(ax=ax, edgecolor=(0, 0, 0, 1), facecolor=(0, 0, 0, 0), linewidths=0.5)
print(grid.head())

# read order file
order = pd.read_csv('E:/DataSet/taxi_data.csv',nrows=100)
print(order.head())

# print()
plt.show()
